Thermodynamics of a 4-site Hubbard model by analytical diagonalization 
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By use of the conservation laws a four-site Hubbard model coupled to a particle bath within 
an external magnetic field in z-direction was diagonalized. The analytical dependence of both the 
eigenvalues and the eigenstates on the interaction strength, the chemical potential and magnetic 
field was calculated. It is demonstrated that the low temperature behaviour is determined by a 
delicate interplay between many-particle states differing in electron number and spin if the electron 
density is away from half-filling. The grand partition sum is calculated and the specific heat, the 
suszeptibility as well as various correlation functions and spectral functions are given in dependence 
of the interaction strength, the electron occupation and the applied magnetic field. Furthermore, for 
both the grand canonical and the canonical ensemble the so called crossin&pohits of the suszeptibility 
are calculated. They confirm the universal value predicted by Vollhardtcl. 
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I. INTRODUCTION 

The Hubbard modefli with the restraint of electron 
hopping to nearest neighbours only is certainly one of the 
simplest many particle models. It was utilized for the de- 
scription of correlated electrons in the intermediate and 
strong coupling regimes. First it was riJ&ed in chemistry 
to describe unsaturated hydrocarbonsoO. In solid state 
physics it was introduced to describe the behaviour of 
narrow-band transition metal compounds, with the main 
emphasis igi-ygn to the explanation of the metal- insulator 
transitiora§0 and the itinerant ferromagnetismH. Later 
on it played a central role in the heavy fermion context, 
and especially the discovery of high T c superconductiv- 
ity in the cuprates renewed the interest, since its involved 
spin dynamics is widely considered as a good candidate 
to understand at least qualitatively the pairing mecha- 
nism in the copper oxide planes, being common to all 
high T c superconductors. It is clear that after about 
four decades of research there exist a vast literature and 
meanwhile also some gopcLreviews dealing with differ- 
ent aspects of the modelETlill, where good reasons to deal 
with that model are given, so I refer to them. In contrast 
to the large amount of work spent to that model, there is 
astonishing few exact knowledge. The initial hope that 
by reducing the description of the correlated electrons to 
a one-band model with hopping restricted to next near- 
est neighbours and correlation restricted to the on-site 
term a model will result, which can be solved analyti- 
cally, has not-yet been fulfilled in general. Nevertheless, 
Lieb and WuEHl calculated the ground-state-energy of the 
one dimensional model with a filling of one electron per 
lattice site by help of the Bethe ansatz. Starting from 
this paper there were later on attempts to get the com- 
plete spectrum and thus the thermodynamics for the half 
filled model and during the last decade also for arbitrary 
electron filling, at least approximately. The state of art 
for the one dimensional model was recently summarized 



by Deguchi et. al.EH. All theories starting from the Bethe 
ansatz result in a set of coupled nonlinear equations, the 
Lieb-Wu equations, which have to be solved numerically 
for every parameter set. Thus all studies of physical 
quantities like susceptibilities or excitation spectra rest- 
ing upon the solution of the Lieb-Wu equations have to be 
done numerically, which often limits the application. The 
transfer matrix method introduced by Shastryta supplies 
an alternate solution, which is in principle exact. The 
drawback for calculating physical quantities is that the 
lowest eigenvalue of the transfer matrix has to be calcu- 
lated, what requires again a considerable numerical ef- 
fort. For the cases where results exist, they .were shown 
to be equivalent to the Bethe ansatz methodic. Another 
way to achieve exact information is the numerical treat- 
ment of finite systems. The techniques by far most itu 
tensively used for this task-are the Lanczos algorithmlij 
and Monte Carlo methodaij. Lanczos calculations. Al- 
beit giving essentially exact results for ground-state and 
low lying eigenvectors and single particle excitation spec- 
tra both methods suffer from the relative small cluster 
size treatable with the contemporary computer facilities, 
since the basis set of vectors grows exponentially with the 
system size. The main drawback in the MC calculations 
for fermions are the relatively large uncertainties due to 
the so called sign problem, although there exist a enor- 
mous number of derivate methods developed to overcome 
this problem. Direct exact diagonalization, e.g. with the 
householder method, is restricted to much smaller prob- 
lems and exact is then always meant numerically. Along 
this way it is in principle impossible to get closed analyt- 
ical expressions in dependence of the model parameters. 
In most cases a numerically exact calculation may be 
enough, but there are problems which need the explicit 
analytical dependence of the eigenvalues and the coeffi- 
cients of the eigenvectors with respect to a suitable basis 
set. The problem of determining the crossing tempera- 
ture of the high temperature specific heat, a problem we 
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address in the third section, is such a case, since it in- 
vokes the third derivative of the thermodynamical poten- 
tial. Another problem, where the analytical knowledge 
of the spectrum is inevitable, is the term-by-term com- 
parison of series expansions with exact solutions, and, 
last but not least, the study of the influence of a mul- 
titude of external parameters, (e.g. in the model under 
consideration the ratio t/U between the hopping energy 
and the on-site correlation energy, the electron density 
n, the temperature, and the external applied fields) is 
considerably simplified, if the eigensystem of an Hamil- 
tonian is known analytically, especially, if the calcula- 
tion aims at the grand canonical partition sum, where 
the chemical potential has to be determined from the 
particle number. To get the explicit dependencies on 
the model parameters it is necessary to diagonalize the 
Hamiltonian in an analytical manner, what restricts the 
calculation to very small systems of course. We have 
chosen a four-point square fo*-jthis task for the following 
reasons: (i) It is the earliestO~c3 and most often stud- 
ied non-trivial object in numerical diagonalization. (ii) 
There were some attempts in the nineties which remained 
incomplete, since they discussed the half filled case only, 
where the eigenvalues were given. The authors speculate 
about hidden symmetries, since they were not able to re- 
solve all degeneraciesEjo. (iii) It is a special case of the 
N-site rings, which are subject of the above cited Bethe 
ansatz and quantum transfer matrix methods, and might 
be useful as a proving tool, (iv) In the following sense the 
four point Hubbard square may serve as a crude model 
to the square lattice: If you neglect every second hopping 
in x- and y-direction respectively you will get a system 
of N/4 independent 4-site-clusters being the topic of this 
paper, see Fig. [I]. The model is thus a little bit bet- 
ter than the so-called atomic limit, were every hopping 
is neglected. The main fault of this cluster-gas model 
is that there is no direct electron transfer between the 
clusters. This is in some sense corrected by applying the 
grand canonical distribution, since it reflects the situa- 
tion where the subsystems are in contact with a particle 
reservoir allowing for particle number fluctuations, thus 
establishing an indirect particle transfer via the bath. 
Very recently the spectral function of the square-lattice 
Hubbard model was calculated by a cluster-perturbation 
method, which used the square-clusters as unperturhed 
system and the inter-square-hopping as perturbatiorc3, 
albeit they diagonalized the cluster numerically, (iv) The 
rigorous solution of any strong interacting many particle 
system is a value independent of its applicability to re- 
alistic systems, at least for pedagogical reasons, (v) The 
limited resources of the available computers, we used for 
symbolic formula handling. 

The paper is organized as follows. In the next two 
sections the diagonalization procedure is explained and 
the many-particle spectrum is discussed. In the fourth 
section we give the complete thermodynamics, based on 
the grand canonical potential. The main focus lies on 
the dependence of the physical quantities on the elec- 



tron density, whereas the comparison of our results to 
the literature for the half-filled case serves as indepen- 
dent proof. In the closing section we discuss, what we 
could learn from our analytical solution. 



II. THE ANALYTICAL DIAGONALIZATION 
PROCEDURE 

We are interested in the eigenvalue sequence for the 
operator 
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Here c^. and Cj CT are the creation and destruction op- 
erators in Wannier representation. The chemical poten- 
tial [i and the magnetic field h in z-direction, which we 
choose to be normal to the square, are introduced to take 
into account the effects of doping and applying external 
magnetic fields. The related spectrum we will call grand- 
canonical, due to its relation to the weights within the 
grand canonical distribution. The lowest level of this 
operator determines the state with the highest weight, 
becoming one for T — > 0, consequently, we will call it the 
grand canonical groundstate, which of course may be dif- 
ferent for different values of n and h. The many particle 
Hilbert space is the direct product of the four single-site 
Hilbert spaces, each containing four states, i.e. the empty 
state, the two single occupied states with spin up and spin 
down, and the double occupied state. Thus we have 256 
states. For example, \u, 0, 2, d) indicates the state which 
has a spin- up electron at the first point, no electron at 
the second position, two electrons at the third, and one 
spin-down electron at the last position. Of course it is 
impossible to diagonalize the resulting 256x256 matrix in 
an analytic manner straight forwardly, instead one has to 
use all the known symmetries or the related conservation 
laws to reduce the matrix to block diagonal form with all 
blocks being at least 4x4. 

The Hubbard model is very good suited for the task 
of analytical diagonalization by help of the conservations 
laws, since it contains several conserved quantities which 
are easy to handle. First of all, there are the two com- 
muting SUI2] algebras of the total spin and the total 
pseudospinEj providing us with four operators commut- 
ing with each other and with the hamiltonian, i.e. 
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We mention here that the total pseudospin is conserved 
only if the hopping matrix elements fulfill the condition 
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what is the case if one restricts e.g. the hopping to next 
neighbours as we supposed in eq. ([I]). The expectation 
value of R 2 is nothing but the difference of the model 
from the half filled case, or in the context of high-T c su- 
perconductivity it is the doping of the model. It can be 
controled by adjusting the chemical potential [i of course. 
Dealing with a special cluster one can use the symmetries 
of it. In our case the actual symmetry is G^ v . Indeed 
there are two works on the four-site square, however re- 
stricted to the half filled case, which tried to use these 
symmetries. Within the firstcj the authors start from a 
parametric expression of the Hubbard interaction, which 
they took obviously from Ref. |27]. But instead of using 
both spin and pseudospin symmetry they used the spin 
symmetry first and then the point symmetry of the four 
site cluster. This way the biggest block in the hamilto- 
nian remains 7x7, what in general prevents an analytical 
solution. Nevertheless they found the roots. More ele- 
gantly the same task was done in Ref. where both 
spin and pseudospin were used first and the translation 
invariance afterwards. This is enough to bring the hamil- 
tonian matrix being 70x70 in the case of the half filled 
model into the wanted block form with none of the re- 
maining blocks being bigger than 4. In both papers the 
authors mentioned a degeneration of the energy levels, 
which they ascribed to some hidden symmetries, a point 
which will be clarified within the present work. 
In this paper we use the same method for arbitrary filling. 
Due to the commutation of R 2 and S 2 with the hamil- 
tonian, we started with the basis of the Hilbert subspace 
with a given value of R 2 or, equivalently, with a given 
number of electrons, and a given value of S 2 . Afterwards 
we calculated the matrix of R 2 and diagonalizcd it result- 
ing in a common eigensystem of R 2 ,S 2 ,R 2 . The next 
step was the diagonalization of S 2 in every subspace for 
given eigenvalues R 2 , S 2 , R 2 thus getting common eigen- 
vectors to R 2 , S 2 , R 2 , S 2 . The last step, if necessary, was 
the calculation of the matrix of the translation operator, 
which fulfills U 4 |/i) = \fi) for a four-site ring and there- 
fore has the eigenvalues 1, i, -1, and -i. We took the 
operator U in the form given in Ref. [2q. 
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are operator representations 
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The operators U, U 2 , U 3 

of the group elements C4 2 , C2 2 , and Cj z L , which form 
together with the 1 standing for E an Abelian subgroup. 
The relation to the conserved momentum is obvious. The 
operator Uy ff permutes two electrons and is therefore 
well-suited to construct the missing group elements. We 
find for 
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Since these operators neither commute with each other 
nor with the U-operators we have to construct new oper- 
ators commuting with all the other operators. The result 
is 
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P(IC 2x ) - P(IC 2y ) 
+P(IC 2a ) - P{IC 2b ) 
P(IC 2x ) - P(IC 2y ) 
-P(IC 2a ) + P(IC 2b ) 
-P(IC 2x ) + P(IC 2y ) 
+P(IC 2a ) - P(IC 2b ) 



-P(IC 2x ) 
-P(IC 2a ) 



P(IC 2y ) 
P{IC 2b ) 



(19) 
(20) 
(21) 
(22) 



GrossetJ constructed further conserved quantities for the 
one-dimensional system with periodic boundary condi- 
tions, which arc in general highly nontrivial due to their 
dependence on the interaction strength. The operator 
given in the first theorem of Ref. 29 reads 



K 1 = Kj + 

K = 2 ( C t+2a C ia - C t C i+2a) 

X ( c i,-cr C i,-0 + C j+1 Ci+^-o-) 



(23) 
(24) 

(25) 



For our four-site ring Kj is zero. In Ref. |9| there an- 
other operator K 2 is given, which does not give new 
results due to our short periodicity length. As men- 
tioned above, for our small system it is not necessary 
to use K 1 , but, since it commutes with all the other 
operators, it must have the same eigensystem like the 
hamiltonian. Since it is an antihermitian operator the 
eigenvalues have to be zero or pure imaginary. Due to its 
rather complicated character we used it among others as 
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a powerful tool for checking the eigensystem. All the cal- 
culations were done by help of algebraic programming. 
The resulting eigenvectors and the related eigenvalues 
were proved by help of the eigenvalue equation. Also 
the matrices for the special case N e — 4 are compared 
to Ref. |24L— where they are given. The matrices Hq,o,i 
and Hi^o^BJ differ in the sign for two off-diagonal ma- 
trix elements, but the eigenvalues are the same and the 
matrices Ho,o,-i an d Ho^.i are interchanged. Since the 
resulting spectrum for the half filled case is not altered 
by these differences, the conclusions drawn in Ref. |24| 
remain true. The complete set of the eigenvalues and 

eigenvectors \m r , r{r + 1), m s , s(s + 1), u, kn,p ++ ,p^ ), 

with m r , m g , s, r, u, kn,p++,p-\ indicating the eigenval- 
ues of the operator-set R 2 , R, S z , S 2 , U, Kn, P ++ , P H , 
in its analytical dependence on the model parameters are 
given in the appendix. Inspecting the spectrum carefully 
we found all the states classified by the eigenvalues of the 
conserved quantities in a unique manner. Thus we clari- 
fied the reason of the mysterious degeneracies mentioned 
in Ref. ||^. Also in Ref. |5J hidden degeneracies remained, 
although not mentioned by the authors. For instance 
the two degenerate states No. 95 and 96 have identical 
eigenvalues of the conserved quantities R z , R, S 2 , S 2 , U, 
which was the set exploited in Ref. ^ij. This degeneracy 
is due toj-the conserved Grosse operator Kn. Heilmann 
and Lief£3 showed for the benzene that there are a mul- 
titude of level crossings. Since they have constructed 
all correlation-independent symmetries they ascribe this 
fact to the existence of unknown U-depend&nt symme- 
tries. We believe, that the Grosse operatorscj are candi- 
dates for the missing symmetry, since they were left out 
in Ref. |2(| For instance, it is obvious that for a six-site 
ring the operator K 1 is U-dependent. Nevertheless, in 
our special case it is not, since the operator Kq vanishes. 
Thus within the four-site model K} represents a further 
U-independent symmetry. By this way we remark that 
there is no contradiction to the non-crossing rule, since 
all the eigenstates differ in the symmetry quantum num- 
bers of mutual commuting operators, a question raised 
in Ref. || 



III. THE GRAND CANONICAL EIGENSYSTEM 

A. The low lying states 

In the following discussion of the spectrum we always 
mean the grand canonical levels, which are numbered for 
referencing according to the appendix, exceptions are ex- 
plicitly indicated. The well-known particle-hole symme- 
try allows to fix the electron number to four if we choose 
/i=U/2. Looking at the lower part of the spectrum of the 
half filled model, as plotted in Fig. |], we find the state 
No. Ill being the groundstate. As already mentioned 
in Ref. |24| this is both a spin and a pseudospin singlet 
(m,.=0 and m s =0) in accordance with the statements of 



Liebls positive-U theoremE3 and the lemma of Shen and 
QiuEl. In Fig. I we show the results for different tem- 
peratures. As the first exited states we find a spin triplet 
(states No. 102, 128, 154), regarding the pseudospin 
it remains a singlet. The next group consists again of 
a singlet-triplet type sequence, where the singlet differs 
from the groundstate by its eigenvalue of the translation 
operator only (it is 1 in comparison with U=-l of the 
groundstate), the triplet states have pure imaginary U 
eigenvalues and are therefore twice degenerated, in agree- 
ment with a statement given in Ref. |24|. It is interesting, 
that the next exited states deviate in the charge. Be- 
ing both spin and pseudospin doublets they are eightfold 
degenerated, if we do not apply neither an external mag- 
netic field nor a deviation of the chemical potential from 
U/2. Since we calculated Fig. || with a small magnetic 
field, the spin degeneracy is lifted. We also find a pseu- 
dospin triplett (not shown in Fig. |^) among the low lying 
states, i.e. the states with No. 22 (n=2), No. 131 (n=4) 
and No. 232 (n=6). A small deviation of /j, from U/2 
will lift the degeneracy of these states, since /1 couples 
to R z in the same manner as the magnetic field couples 
to S z . But in between the pseudospin singlet and triplet 
states we find a further pseudospin-singlet-spin-quintet- 
structure and pseudospin-doublet-spin-quartet structure. 
It was quite surprising for us, that the low lying spec- 
trum of the half-filled model contains states with electron 
numbers differing from four by one and two. For finite 
temperatures and deviation from half-filling we have to 
determine the chemical potential from the N(fi) curves. 
Due to the analytically given expression, this is in prin- 
ciple possible in any accuracy. For example, the chemi- 
cal potential used in the Fig. [|, where we depicted the 
lowest part of the spectrum for (N e )=3, results in a de- 
viation of the electron number from 3 which is less than 
10~ 12 for T = 1CP 7 <. If the temperature reaches zero, 
we have a step curve. Surprisingly the model steps from 
N e =2 to N e =4 with no plateau at N e =3. This corre- 
sponds to the interesting fact that there is an extreme 
small /i-window, where a n=3 state is groundstate. This 
is demonstrated by help of the inset of Fig. ^ where we 
plotted the groundstate in dependence of the chemical 
potential. Thus it is interesting to look at the spectrum 
with N e =3 also. Since we work with the grand-canonical 
potential we have to determine first the chemical poten- 
tial. Fig. |5| makes clear that we find again states with 
different electron numbers among the lowest levels. The 
degeneration of the pseudospin-triplet is now lifted due to 
the large deviation of the chemical potential from U/2, 
and its lowest state (No. 22) is the first exited state 
above the groundstate, which is No. 46 and No. 50. 
Both states belong to the N e =3 pseudospin doublets in 
the half-filled model. As mentioned above the degeneracy 
of the groundstate is due to the imaginary eigenvalue of 
the translation operator. Among the low lying states the 
No. Ill, being the groundstate of the half-filled case, 
remained. If we do not apply a small magnetic field, 
the states No. 46 and No. 50 are degenerated with the 
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states No. 70 and 74, and so do the states No. 22 and 
No. Ill, which become the groundstates of the model 
for lower correlation, although their occupation eigen- 
value is not 3. This situation is depicted in Fig. @. Due 
to the interplay of the four s=l/2 states, belonging to 
N e =3, and the two s=0 states belonging to N e =2 and 
N e =4 respectively one has to expect a complex magnetic 
behaviour of the model. For example a small magnetic 
field lifts the degeneracy of the N e =3 states whereby the 
states No 46 and 50 are lowered below the states No. 
22 and No. 111. This process is shown in Fig. |7|. By 
further increasing the correlation energy U we found a 
spin change of the groundstate, since the S=3/2 states 
no. 38, 58, 82, and 90 become the grcmndstates. This 
was first discovered by Callaway et.al.E3 via numerical 
diagonalization. We add that the groundstate is fourfold 
degenerated, and the exact U value, where the transition 
happens is U/t=4(2+v / 7), what differs slightly from the 
numerical value 1/0.055 given in Ref. The difference 
is in the second digit and may serve as an estimation of 
the accuracy of their calculation. For N e =5 they did not 
find such a change in the groundstate spin, what is con- 
tradictory to the electron-hole symmetry, which was oth- 
erwise explicitly mentioned in their paper - and of course: 
the same transition happens at the same U value. 



B. The contents of double occupied states 

A further question of interest is whether projecting out 
the double-occupied states from beginning is a good ap- 
proximation. Since we have all the eigenstates in analyt- 
ical form, it was easy to calculate the contents of double 
occupied states in each state. The U-dependence of the 
admixing of double occupied states to the lowest lying 
states is depicted in Fig. g. Among the considered states 
the groundstate has the highest content for bigger U, be- 
ing about ten per cent for U=10t and falling to about 
two per cent for U=30t. If the correlation energy U is 
lower than the bandwidth the double occupied states ad- 
mix considerably. Thus one has to be careful with the 
projection, especially if one aims at a description of the 
Hubbard model in the region were correlation energy and 
bandwidth are of the same order. 



with Z being the grand canonical partition sum, and E m 
being the grand canonical levels. The most interesting 
function due to its relation to spectroscopy we get with 
B = Ciu and A = c+.. In Fig. g we show the results in 
dependence of U for the case of half-filling. The most oc- 
ular feature of the spectral function is the opening of the 
gap of width U with increasing correlation strength. This 
is in qualitative accordance with the common picture of 
the Hubbard-Mott transition. 



IV. THERMODYNAMICS 

Once calculated the spectrum in an analytical way it 
is straight forward to get the grand canonical potential 
and the complete thermodynamics. For the four point 
cluster there exist several papers, which got the ther- 
modynamics by numerical diagonalization. Heinig apd 
Monecke calculated for this model the specifp, heatEH, 
the susceptibility and local magnetic momentsta for the 
half filled case. Shiba and Pincus did the same within a 
comprehensive work devotecLto the one dimensional half- 
filled-band Hubbard modelHl. The discrepancies of these 
two works are discussed in a paper of Cabib and Kaplan, 
who calculated in contrast, both the canonical and the 
grand canonical ensembloEj. They especially discussed 
the specific heat and the spin-spin correlation functions, 
where the main focus was again the half-filled case. In 
the following we give our results for both the half-filled 
case (for comparison) and the doped case. 



A. Thermodynamic density of states 

In their paper@ Cabib and Kaplan argued that the 
grand canonical ensemble for a finite N system reflects 
the infinite N case much better than the canonical ensem- 
ble, especially for the large U case. We agree with that 
point of view, especially regarding the thermodynamical 
density of states (TDOS). We calculated this quantity, 
i.e. the derivation of the particle number with respect to 
the chemical potential. 



C. The spectral function 



D T (n) = 



dNjjii) 



(27) 



By help of the eigenstates and the spectrum it is possi- 
ble to calculate the spectral function for arbitrary opera- 
tors. Since the <5-like peaks result in discontinuous curves 
being inconvenient for the graphics we substituted the 5- 
functions by Lorentz distributions, i.e. we defined 
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7T S 2 + (to - E m + E n f 



(26) 



In Fig. [L0| we depicted the results for different values 
of U at the relatively high temperature of T = O.lt for 
graphical reasons. It is obvious, that the TDOS changes 
its character qualitatively with increasing U starting 
from a three-peak structure for small on site correlation 
changing to an eight-peak structure at large correlation. 
Around U = It there opens a gap around /i = U/2. 
This means, adding a further electron to the half-filled 
model needs an energy of order U. We interpret this, in 
the same sense as above for the spectral function, as the 



5 



variant of the "Hubbard-Mott-transition" for our small 
model. Of course, the "critical U", where the gap opens, 
depends on the temperature. For T = the gap opens 
for infinitesimal small values of the correlation. This is a 
special case for the same statement, proved for the one- 
dimensional Hubbard model by help of the [exact solution 
for the groundstate given by Lieb and WutHl. 



B. Specific heat 



The temperature dependence was controversial in 
Refs. [l7| , |l9| , po| . For the half-filled case we find exactly 
the same features as given in Ref. |2(]. Cabib et. al. 
stated the qualitative picture given in Ref. [l?] to be right 
and that of Ref. |l9| to be wrong for the case U<6. Un- 
fortunately a quantitative comparison with Ref. [l?] was 
impossible due to inconsistencies in that paper. In Fig. 
|IT| we show the specific heat curves in dependence on the 
temperature. The curves with U=0.7t and U=8t respec- 
tively are exactly the same as given in Fig. 1 of Ref. |2(] 
and thus confirms this paper and serves as independent 
proof of our calculation. Furthermore, the change from 
the double peak shape to the triple peak structure is due 
to the behaviour of state 119, what is evident from Fig. 
|l2| . Variation of the electron density causes dramatic 
changes of the specific heat. In Fig. [l3| we calculated for 
a constant electron density and constant magnetic field 
(for simplicity set to zero). The shape of the specific heat 
curves is qualitatively affected if the electron density is 
changed. In the right column of Fig. [Ds] the ratio of cor- 
relation energy U and band width is chosen to 1. This 
is the most interesting case, since it is not accessible by 
means of perturbation theory. We see that with devia- 
tion from half-filling the small peak at ke=0.1t vanishes, 
whereas a low energy peak develops at ksT=0.002t, be- 
ing maximum at N e =3, i.e. n e =0.75. The reason is the 
small distance between the states 46, 50, 70, 74 belonging 
to N e =3, and the state 111, belonging to N e =4, and state 
22, belonging to N e =2, as may be seen from Fig. |5[ Since 
the N e =3 groundstate is magnetic and the groundstates 
for N e =2 and N e =4 respectively are not, these excita- 
tions are both spin excitations and charge excitations. 
In the right column of Fig. |l3| the correlation strength 
is twice the bandwidth. In that case we find no small 
energy peak for N e deviating from 4. Instead the double 
peak structure smears out due to the increasing influ- 
ence of magnetic excitations which are maximum if the 
electron number is fixed to three. Further decreasing the 
electron number results in vanishing of the double peak 
structure. The reason is that the number of accessible 
states with double occupation is reduced drastically. 



C. The crossing points 

A further subject of interest are the so called crossing 
points, which attracted much interest recently, since they 
are present in several strong correlated systems. Voll- 
hardt clarified this pointEl. In a second paperL^I he and 
the coworkers specified to Hubbard models and found 
a (nearly) universal crossing point of the C(T,U)/ks- 
curves. Since the arguments in Ref. ^ never use the 
thermodynamic limit, we feel that they should work for 
our small model also. In Fig. [l4| we depicted the spe- 
cific heat curves nearly in the same fashion as was done 
in Fig. 1 (b) of Ref. |33|. The similarity is ocularly. 
The most striking point the authors made in Ref. |33| 
is the discovery of a nearly universal value of about 
0.34 of the specific heat at the high temperature cross- 
ing point for small U. We find for our model a value of 
the crossing temperature T*(U^0)/t=1.334820204t and 
C(U^0,T^)=0.3351877115 Nk B . All the Hubbard sys- 
tems considered in Ref. [53] are in the thermodynamic 
limit at half-filling. Whether our small cluster is com- 
parable to infinite systems is a difficult question. Usu- 
ally the physical quantities calculated from the canonical 
and grand-canonical potentials differ in small systems. 
Therefore this difference may serve as a rough indicator 
how relevant the results are for the thermodynamic limit. 
In Fig. [[5] we plotted the specific heat for the canonical 
ensemble. We found again a high temperature crossing 
point, whereas there is no sharp crossing point around 
k B T/t=0.8, in contrast to Fig. |l|. The related spe- 
cific heat value is C(U^0,T*)=0.3744085787 Nk s with 
T*=0.940347332t. This is only a little bit bigger than the 
grand-canonical value. Nevertheless the big difference in 
T* shows that there are differences at all. The crossing 
points at different U values were determined according 
to Ref. [l] by help of the equation 



dC{T,U,y) 
dU 



= with fj, = U/2 



(28) 



(29) 



for the grand canonical ensemble and by 
dC{T,U) 

du 

for the canonical ensemble with N e =4 respectively. These 
functions significantly differ from each other, as one can 
see from Fig. [l6|. 



D. Susceptibility 

The isothermal susceptibility shows the interplay be- 
tween the low-lying nonmagnetic states and the low-lying 
magnetic states. As shown above, the fact whether the 
groundstate is magnetic or not is dependent on the oc- 
cupation. That is why we find an increasing peak due 
to the low-energy magnetic excitations if we decrease the 
electron number from four to three, which vanishes, if we 
subsequently decrease the number of electrons towards 
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two. This is illustrated in Fig. [l7| In the strong cor- 
related case (U/t=8) we find a complete different be- 
haviour. Now the distance between the states with differ- 
ent electron occupation is large. Consequently, we do not 
find any low-energy peaks as visible in Fig. Instead 
we see a divergent behaviour due to the S z degenerated 
groundstate if the electron number deviates from four. 
It is well believed, that the Hubbard model in three di- 
mensions exhibits antiferromagnetism near the half-filled 
limit and for low temperatures. In our small model long 
range order does not exist, nevertheless we can ask how 
the correlation functions depend upon change of the cor- 
relation strength, of temperature and electaon filling. We 
adopt the notation of Cabib and KaplanEJ for the spin- 
spin correlation functions 

L n = 4(Sf Sf +n ) (30) 

If not otherwise mentioned we always take the expecta- 
tion with respect to the grand canonical ensemble. By 
help of the functions L n one can calculate the suscepti- 
bility for our small cluster according to 

X T = 4/3(L + 2ii + L 2 ) - 4/3(S 2 ) . (31) 

The susceptibility calculated from eq. |l| has to agree 
with the expression we got from the differentiation of the 
thermodynamical potential. We used this fact to check 
the calculated correlation functions. In Fig. [l^ the three 
different spin-spin-correlation functions are shown in de- 
pendence on the temperature for the half-filled case. The 
model shows the well known antiferromagnetic coupling 
of adjacent spins. Of course this breaks down if the tem- 
perature is increased. The breakdown shifts to lower tem- 
peratures when the strong coupling regime is reached. 
This is exactly what is expected, due to the fact that for 
large values of U the model can be mapped onto an an- 
tiferromagnetic Heisenberg model with the coupling con- 
stant being cx 1/U. Obviously the U=4t case and the 
U=8t case have nearly the same breakdown tempera- 
ture, fortifying the above stated fact, that U=8t, i.e. U 
is twice the unperturbed bandwidth, belongs not to the 
strong coupling limit. It is exactly this parameter regime 
where we find an increase of the local magnetic moment 
with increasing temperature before it breaks down to its 
high temperature limit. The strong correlation behaviour 
of the model is shown in Fig. |2fj. The most striking 
fact is the change from antiferromagnetic correlation to 
ferromagnetic correlation if the electron density is low- 
ered from one to three quarters. Also the next nearest 
neighbour correlation changes its sign when the electron 
occupation changes from three towards two. While for 
the neutral model the nearest neighbours are antiferro- 
magnetically correlated, the next nearest neighbours are 
ferromagnetic. For three electrons both correlation func- 
tions are ferromagnetic and for two electrons they are 
both antiferromagnetic. The temperature dependence of 
the density-density correlations is depicted for the neu- 
tral model in Fig. ||l]. Regarding (nini) = (ni + 2d) 



with d being the double occupancy operator d = n^n^ it 
is obvious how the onsite electron correlation suppresses 
the probability to find two electrons at the same lattice 
site as long, as the correlation energy is small with respect 
to the temperature. Especially the correlation function 
(111112) makes evident that the tendency to charge or- 
dering is lowered. We remark, that the TDOS may be 
calculated via the density-density correlation functions, 
giving one more possibility to prove the inner consistency 
of the derived formulae. In Fig. ^] we depicted the 
nearest-neighbour density correlation for strong correla- 
tion in dependence of the electron filling of the model. 
With decreasing electron number the electrons get more 
space. Due to the repulsive electron-electron correlation 
the electrons tend to increase their distance. Therefore 
the nearest-neighbour density-density correlation is re- 
duced considerably. 



V. DISCUSSION 

The benefit one gets from an analytical solution of a 
model is the easy access to its complete physics. Due to 
the limited space, we had to restrict the presentation of 
the results to a few points. We focused here mainly to 
two points. The first one was the study of eigenvalues 
and related eigenstates in dependence on the correlation 
strength. This way we clarified the open questions about 
possibly hidden symmetries. It was shown, that every 
state is distinguished from each other by using the two 
SU[2] symmetries of spin and pseudospin, the translation 
symmetry, and the eigenvalues of the Grosse operator 
Ki. Since therefore all the eigenstates are off different 
symmetries, there remains no problem with the multi- 
tude of level crossings present within the model. The low- 
est states in the half-filled model yield the singlet-triplet 
picture, widely used in the context of high-T c supercon- 
ductivity. In the half-filled model the charge excitations 
are a little bit higher in energy, but, these excitations 
are lowered considerably in energy, if the electron num- 
ber is driven towards three by appropriate tuning of the 
chemical potential, thus enlightening the complex inter- 
play of spin and charge (more exactly: pseudospin) de- 
grees of freedom. We were really surprised that the states 
with adjacent electron numbers are still lower in energy 
if three is chosen for the mean occupation, even if the 
correlation equals the bandwidth. This situation is com- 
pletely changed if one applies a (relatively) small mag- 
netic field or if the ratio U/t is increased slightly ( what 
may be done by help of pressure in real systems). Since 
our model is a special case of the n-site Hubbard rings 
which are the pet of the Bethe ansatzers, there must be 
a correspondence between the Bethe ansatz eigenstates 
and the eigenstates given in the Appendix. However, it 
is not obvious, how the quantum numbers characteriz- 
ing the states in our solution correspond to the quantum 
numbers of the theories based on the Bethe ansatz. Two 
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of them, the z-component of spin and pseudospin respec- 
tively are used within the Bethe ansatz, but regarding 
the other quantum numbers of the Bethe ansatz, i.e the 
charge momenta and the spin rapidities, we have no idea 
how they are related to the remaining quantum numbers 
of our model, i.e. the modules of spin and pseudospin, 
the eigenvalue of the translation operator, and the eigen- 
value of the Grosse operator. The complicated numerical 
solution of the Lieb-Wu equations delivering the charge 
momenta and spin rapidities prevents a deeper insight, 
so it remains an open question. The interesting problem 
whether the eigenstates may be classified as k strings 
and k-A strings as Deguchi et. al.Ej did with respect to 
solutions of the Lieb-Wu_equations following the string 
hypothesis of TakahashiEll, or in the way as was done in 
a very recent papercj remains a unsolved question for 
the time being. A second point we were interested in, is 
whether the model gives some hints towards the metal- 
insulator transition. Both the spectral density and the 
thermodynamical density of states show the expected be- 
haviour, i.e. the opening of a gap of width U for T>0, 
so there is nothing spectacular in this direction. The 
general theorems for the Hubbard model are fulfilled of 
course. Our analytical solution allowed a detailed study 
of the crossing points of the specific heat. Our results 
confirm the explanation given by Vollhardt et. al. com- 
pletely At the first glance this seems astonishing, since 
we have a level spectrum instead of the various contin- 
uous DOS' used in Ref. [33]. On the other hand, if one 
deals with the high-temperature crossing point the tem- 
perature smoothes the level structure in some sense. Fur- 
thermore, the argument, given in Ref. |33|, that the energy 
scale for T is essentially t if we are at higher tempera- 
tures and of the order of the singlet-triplet excitation at 
lower temperatures holds for our model also. The com- 
parison of the canonical and grand canonical ensemble 
at half filling exhibit differences regarding the tempera- 
ture dependence of the derivatives of specific heat with 
respect to the correlation strength. The reason is that a 
"sampling" of the DOS is much poorer with 72 states in- 
stead of 256. Thus, it is not unexpected that the values 
of C + differ slightly from the value calculated in Ref. |33| 
for an infinite linear chain by help of second order per- 
turbation theory in U/t. Considering the specific heat 
for deviations of the mean occupation number from half- 
filling we find another much lower energy scale, due to 
the fact that states with different electron numbers are 
nearly degenerated if we fix the mean electron number 
to three. Regarding the magnetic behaviour we see the 
interplay between the even-occupied clusters, which are 
nonmagnetic and exhibit low-lying singlet-triplet excita- 
tions and the odd occupied clusters which exhibit a mag- 
netic degenerated groundstate with doublet-quartet ex- 
citations if we have N e =3 or N e =5. By fixing the mean 
occupation number to non-integers we enforce a mixture 
of both cluster groups yielding a diverging susceptibil- 
ity due to the contents of magnetic clusters. The spin- 
spin correlations confirm the common accepted picture, 



that the strong on-site correlation favours antiferromag- 
netism instead ferromagnetism. Otherwise, next-nearest 
neighbours are coupled ferromagnetic. This holds for the 
half filled model. This situation is changed completely 
if the mean occupation number is driven towards three. 
There we find ferromagnetic coupling between adjacent 
spins for low temperature. A completely unexpected be- 
haviour is the change to an antiferromagnetic coupling 
(albeit small) when the temperature is increased and be- 
comes of the order of the hopping matrix element. For 
higher temperatures the off-site spin-spin correlation van- 
ishes. Since we have a small cluster long range mag- 
netic order is impossible, nevertheless, if one regards the 
nearest neighbour correlation as an hint to the magnetic 
phases realized in an infinite system, our findings support 
the common accepted phase diagram for the (simple cu- 
bic) three-dimensional Hubbard model. The suppression 
of electron motion with increasing correlation strength 
is evident from the density-density-correlations. If U is 
large enough we find exact one electron at every site. Mo- 
tion seems forbidden, what hints towards an insulator. 
Regarding the square lattice we believe that it is more 
appropriate, with reference to the conserved quantities 
of the Hubbard model, to speak about spin-pseudospin 
separation instead of spin-charge separation, since the 
"charge" is the "z-component" of the pseudospin only. 
From models in magnetism it is common knowledge that 
quantum fluctuations reduce the magnetic moment, due 
to the fact that the total spin is conserved, while the lo- 
cal spin is usually not. In the Hubbard model the same 
holds for the pseudospin also. If the hopping between the 
clusters is treated by a perturbation theory as was done 
in Ref. |25|, the calculated cluster levels will form bands, 
and the related cluster-spins and -pseudospins will not 
be conserved furthermore, thus "pseudospin- waves" will 
arise as well as spin-waves. The undoped model has a 
vanishing total pseudospin. Reducing the number of elec- 
trons from one per lattice cite creates a certain number 
of empty places. At these sites a local spin is missing 
now and a local pseudospin is created. Further dop- 
ing dilutes the spins and increases the distance, what 
reduces the spin-spin interaction. Otherwise the number 
of pseudospins increases; their mean distance is reduced 
and the pseudospin-pseudospin interaction is increased, 
thus favouring long-range order of the pseudospins. As 
we have shown, among the lowest eigenstates of our clus- 
ter we find states with different spins and pseudospins. 
Since the related states are many-particle states the spins 
and the pseudospins are not coupled to individual elec- 
trons or holes. If these many-particle states become itin- 
erant by switching on the inter-cluster-hopping, there is 
no reason, that the resulting bands must exhibit identical 
shapes, or in other words that the resulting quasiparti- 
cles have the same mass. On the contrary, one would 
expect, that a quasiparticle made from a cluster-state 
with three electrons and spin 1/2 will differ significantly 
from a quasiparticle made from a spinless cluster-state 
with four (or two) electrons, both in the effective mass 
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and in the velocity of propagation. This is exactly, what 
was calculated in Ref. |2f] and what is called spin-charge 
separation. Finally, we mention that the given solution is 
applicable to the four-site Hubbard model with attractive 
interaction as well, simply by applying the Shiba mapta 
— * ai and Cij — * a^j, which exchanges the spin with 
the pseudospin. The tables given in the appendix apply 
as well for negative U hamiltonian. The only task is to 
replace U — > —U, h + U/2 — > fi, and fx + U/2 — > h for the 
energy eigenvalues and to exchange the spin and pseu- 
dospin quantum numbers in the assignment of the eigen- 
vectors. All the presented figures, may be calculated for 
the attractive model as well. There is a multitude of other 
interesting features, e.g. the transport properties within 
the cluster, which are nevertheless easy to calculate by 
help of the analytically given eigenstates and eigenvec- 
tors. We abstain from presenting the related figures to 
prevent overloading of the present paper. In summary 
we may say that the analytical solution of the four-site 
Hubbard model for arbitrary electron filling, arbitrary 
interaction strength and arbitrary temperature is highly 
non-trivial. It allows easy access to all interesting fea- 
tures of the model and may serve as a reference object for 
various numerical or perturbation methods dealing with 
more complex models of strong electron correlation. 



ACKNOWLEDGMENTS 

I would like to thank S.-L. Drechsler and R. Hayn for 
usefull discussion, J. Monecke for some helpful informa- 
tion on their early work, and P. Vollhardt for encouraging 
remarks on the crossing point calculation. 



APPENDIX: EIGENVECTORS AND 
EIGENVALUES 



Here we present some eigenstates and the com- 
plete spectrum. Since by far the most papers deal 
with the groundstatc of the Hubbard model we re- 
strict to presentation of the groundstate of the neu- 
tral model with N e =4, i.e. state 111, and one of the 
groundstates for N e = 3, i.e. state 46. The an- 
alytical form of the other 254 eigenvectors the inter- 
ested reader may find at the web-page www.physik.tu- 
dresden.de/itp/members/schumann/rescarch.html. 

$ m = Ci (|0022) - |0220) - |2002) + |2200>) 

+C" 2 (-|02du) + \02ud) + \0du2) - \0ud2)) 
+C 2 (-|20du) + |20ud) + \2du0) - \2ud0)) 
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In following Tables |j- XIH we give the results for the 
whole spectrum. The eigenvalues of the operators 
R z ,R,S Z ,S 2 ,U,Kn, P++,P+-, H are assigned as 
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m r , m s , r, s, it, kit, p++, Ph , £■ In the first col- 
umn we number the eigenstates. In the second col- 
umn the energy eigenstates are given in the form 
\m r , m s , r(r + 1), s(s + 1), u, ku,p ++ ,p-\ — ). The related 
energy eigenvalues in dependence on the parameters t 
and U are given in the third column. For the grand 
canonical eigenvalues these energy eigenvalues have to 
be completed by —pN e — hm s , therefore we give the re- 
lated values in the subheads of the tables. In the fourth 
column we calculated the numbers for the energy for 
t=l, U=5, fi = U/2, and h = 0.01, i.e. the half filled 
case. The small magnetic field was introduced to lift 
spin-degeneracy. 
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FIG. 1. By replacing every second hopping line by an 
indirect hopping via the bath the square-lattice (A) is trans- 
formed to the square-cluster gas (B). 

FIG. 2. The low lying states of the half-filled model. The 
parameters are t=l, U=5, h=0.01, (iV e )=4. The main quan- 
tum numbers are indicated at the levels. The length indicates 
the degeneracies. 

FIG. 3. The particle number in dependence on the chem- 
ical potential. The parameters are t=l, U=5, h=0.01. 

FIG. 4. The grand canonical groundstate in dependence 
on the chemical potential. The parameters are t=l, U=5, 
h=0.01 . 

FIG. 5. The low lying states for {N e )=3. The parameters 
are t=l, U=5, h=0.01 (h=0 for the dotted curve), k s T=l(T 7 . 

FIG. 6. The lowest states for (iV e }=3 and U=4t and U=5t 
respectively. The remaining parameters are t=l, h=0 

FIG. 7. The lowest states for <iV e )=3 and h=0t and 
h=0.01t respectively. The remaining parameters are t=l, 
U=4t 

FIG. 8. Contents of double-occupied states. The parame- 
ters are t=l, h=0, {N e )=4. The operator P projects out the 
double-occupied states. 

FIG. 9. The spectral function J + . The parameters 
are t=l, h=0, (iV e }=4. 

FIG. 10. The thermodynamical density of states. The 
parameters are t=l, h=0, k_BT=0.1t. 

FIG. 11. The specific heat in dependence on the tempera- 
ture for different values of the correlation strength calculated 
from the grand canonical partition sum. The other parame- 
ters are t=l, h=0, (iV e )=4. 

FIG. 12. The U-dependence of the low-lying levels for the 
half-filled model. 



FIG. 14. The specific heat in dependence on the tempera- 
ture for the grand-canonical ensemble at half-filling. The up- 
per picture shows at a large scale the crossing points, whereas 
the inset shows the high temperature crossing point in fine 
scale. The lower plot demonstrates that there is no definite 
high-temperature crossing point for large U values. The pa- 
rameters are t=l, h=0, and the U values are indicated at the 
curves. 

FIG. 15. The specific heat in dependence on the temper- 
ature for the canonical ensemble and half-filling. The inset 
shows the high-temperature crossing point in fine scale. The 
U values are the same as in fig. |lj. 

FIG. 16. The partial derivative of the specific heat with 
respect to U for the canonical and grand-canonical ensemble 
respectively. 

FIG. 17. The temperature dependence of the isothermal 
susceptibility for t=l and U=4. The parameters at the curves 
indicate the electron number. 

FIG. 18. The temperature dependence of the isothermal 
susceptibility for t=l and U=8. The parameters at the curves 
indicate the electron number. 

FIG. 19. The Spin- Spin- Correlations in dependence on 
the temperature for the correlation parameters t=l, U=l, 
U=4, and U=8. The parameters at the curves indicate the 
electron number. 

FIG. 20. The nearest neighbour spin-spin-correlation 
function (S1S2) for strong correlation U=40 in dependence 
on the temperature. The parameters at the curves indicate 
the electron number. The inset shows the density dependence 
of the spin-spin correlation in the groundstate. 

FIG. 21. The density-density-correlations in dependence 
on the temperature for the half-filled model and for the cor- 
relation parameters t=l, U=l, U=4, and U=8t. 

FIG. 22. The nearest neighbour den- 

sity-density-correlation function (111112} for strong correlation 
U=40 in dependence on the temperature. The parameters at 
the curves indicate the electron number. The inset shows the 
density dependence of the nearest-neighbour density-density 
correlation in the groundstate. 



FIG. 13. The specific heat in dependence on the tempera- 
ture for different electron filling. The number inside gives the 
number of electrons. The parameters are t=l, h=0, and U=4 
for the left column and U=8 for the right column respectively. 
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TABLE I. Eigenkets and eigenvalues for N c =0. 



Number Eigenstate Energy eigenvalue 



N e =0 m s =0 r(r+l)=6 s(s+l)=0 



1 


-2,0,6,0,1,0,4,0) 







TABLE II. 


Eigenkets and eigenvalues for N e =l. 




Number 




Energy eigenvalue Value for U/t=5, /i=U/2, and h=0.01t 




AT 1 

N e =l 


m s = I r(r+l) = -f s(s+l) = § 




2 


j 3 1 15 3 lOO A \ 
\ 2 ' 2 ' ~4~ ' 4 ' ' > > / 


-2t 


-4.51 


3 


' " TT * TT ■ , - "7 - " £ ■ w / - ■ 0/ 

2 2 4 4 77/ 





-2.51 


4 







-2.51 


5 


1-3115 3 1 4 0) 


2t 


-0.51 




N e = l 


m s =-± r(r+l)=f s(s+l) = § 




6 


1 _ 3 _1 15 3 _ x q q _ 4 \ 

1 2 ' 2 ' 4 ' 4 ' ' ' ' / 


-2i 


-4.49 


7 


1 — 2'~2 ,- 4 - '4'~^'^^'^'^) 





-2.49 


8 


|-|,-i,f,f,i,-2i,0,0) 





-2.49 


9 


|_3_1153 1 q 4 q\ 

2 ' 2 ' 4 ' 4 ' ' ' ' / 


2t 


-0.49 
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TABLE III. Eigenkets and eigenvalues for N e =2. The abbreviation a is denned as a — arccos(^-3-^ — %f/[ + \- ) )• 



Number 


Eigenstate 


Energy eigenvalue 


Value for U/t=5, £t=U/2, and h=0.01t 


N c =2 m s = l r(r+l)=2 s(s+l)=2 


10 


1-1,1,2,2,-1,0,0,-4) 





-5.02 


11 


| -1,1, 2,2, -i, -2 1, 0,0) 


-2t 


-7.02 


12 


- l,l,2,2,-i,2t,0,0) 


2t 


-3.02 


13 


- l,l,2,2,t,-2t,0,0> 


2t 


-3.02 


14 


| -1,1,2,2,1,2*, 0,0) 


-2t 


-7.02 


15 


1-1,1,2,2,1,0,-4,0) 





-5.02 


N c =2 m s =0 r(r+l)=2 s(s+l)=0 


16 


- 1,0,2,0, -1,-2 iV2, 0,0) 


-5. 


17 


| - 1,0,2,0, -1,2 i\/2, 0,0) 





-5. 


18 


- l,0,2,0,-i,0,0,0) 


U \/l6t 2 + U 2 
2 2 


-5.70156 


19 
20 
21 


j - l,0,2,0,-t,0,0,0) 
| - 1,0,2,0,1,0,0,0) 
- l,0,2,0,t,0,0,0> 


U | V 16 t 2 + u 2 

2 2 

u V 16 t 2 +u 2 

2 2 

U | \/l6t 2 + U 2 

2 ' 2 


0.701562 
-5.70156 
0.701562 


22 


j - 1,0,2,0, 1,0,4,0) 


f - | V48i 2 + t/ 2 cos(f) 


-8.34789 


23 


1-1,0,2,0,1,0,4,0) 


f + | V48t 2 + (7 2 cos(^p) 


1.5136 


24 


1-1,0,2,0,1,0,4,0) 


§- + § V48t 2 + I7 2 cos(^) 

it 2 u u 3 

a = arccos a ^—r 

( 16 t 2 , £7 2 \ 2 
3 + 9 j 


-3.16571 


N c =2 m s =0 r(r+l)=2 s(s+l)=2 


25 


- 1,0,2,2,-1,0,0,-4) 





-5. 


26 


- 1,0,2,2, -i, -2 i, 0,0) 


-2* 


-7. 


27 


- l,0,2,2,-t,2i,0,0> 


2t 


-3. 


28 


- l,0,2,2,i, -2 i, 0,0) 


2t 


-3. 


29 


- l,0,2,2,i,2z,0,0) 


-2t 


-7. 


30 


j- 1,0, 2, 2, 1,0, -4,0) 





-5. 


N c =2 m s =0 r(r+l)=6 s(s+l)=0 


31 


| - 1,0,6,0,-1,0,0,-4) 


£/ 





N c =2 m s = -l r(r+l)=2 s(s+l)=2 


32 


1-1,-1,2,2,-1,0,0,-4) 





-4.98 


33 


j -1,-1, 2,2, -i,-2i, 0,0) 


-2t 


-6.98 


34 


| - l,-l,2,2,-i,2i,0,0) 


2t 


-2.98 


35 


- l,-l,2,2,i,-2i,0,0) 


2t 


-2.98 


36 


- l,-l,2,2,i,2i,0,0) 


-2t 


-6.98 


37 


1-1,-1,2,2,1,0, -4,0) 





-4.98 
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TABLE IV. Eigenkets and eigenvalues for N c =3, spin-up states 



Number 


Eigenstate 


Energy eigenvalue Value for U/t=5, /tt=U/2, and h=0.01t 


N c =3 m s =§ r(r+l) = | s (s+l) = f 



f ,-1,0,0,4) 
f ,-i,-2t,0,0> 
f ,t,2i,0,0> 



f ,1,0, -4,0) 



-It 



It 



-9.53 
-7.53 
-7.53 
-5.53 



Nc 



=3 m s = ^ r(r+l) = | s(s+l) = f 



|,-l,-t>/3,0,0> 
|,-l,-iV3,0,0> 
§,-MV3,0,0) 
|,-l,i>/3,0,0> 

|,-t,t,0,0> 

|,-t,t,0,0> 



t/ _ 

2 

2 ~ 

t/ _ 
2 

2 ~ 





-44 [7+(7 2 


2 


Viet 2 


-4i C7+(7 2 


2 


Viet 2 


-44 [7+(7 2 


2 


Vwt 2 


-it U+U 2 



|,-t,t,0,0> 

I,-*,*, 0,0) 
I,*,-*, o,o) 



[,t,-»,0,0> 
},»,-», 0,0) 



- + 

2 ^ 



2 



- + 

2 ^ 



f,t,-t,0,0> 

1,1,-1^3,0,0) 

|,1,-»V3,0,0> 

f,l,iV3,0,0) 

1,1,2^3,0,0) 



•J 32 t 2 + U 2 +4 \Jiii t 4 +3 t 2 u 2 
2 

■J 32 t 2 + (7 2 +4 i 4 +3 t 2 (7 2 

2 

U y 32i 2 + (7 2 -4 Vf54t 4 +3t 2 C/ 2 ^ 
~ 2 
[/ y 32i 2 + !7 2 -4 v / 64t 4 +3t 2 U 2 
~ i 2 

y 7 32 t 2 + U 2 +i \Jiii t 4 +3 t 2 u 2 
2 

\/ 32i 2 + (7 2 +4 v / 64t 4 +3t 2 U 2 
2 

^/ 32i 2 + (7 2 -4 v / 64t 4 +3t 2 C/ 2 
2 

32i 2 + (7 2 -4 V64t 4 +3t 2 C/ 2 
2 

Vie t 2 +4t u+u 2 

2 

Vl6t 2 +4tC/+l/ 2 ' 



u _ 

2 

2 + 



u 

2 



2 + 

2 ^ 



Vl6t 2 +4tt/+C/ 2 
2 

Vl6t 2 +4tl/+l/ 2 



-7.30129 
-2.71871 
-7.30129 
-2.71871 

-10.1129 

0.0929233 

-6.57849 

-3.44151 

-10.1129 

0.0929233 

-6.57849 

-3.44151 
-8.91512 
-1.10488 
-8.91512 
-1.10488 



N c =3 m s = i r(r+l) = § s(s+l) = 



f ,-1,0,0,4) 
f ,-i,-2t,0,0> 
f ,t,2»,0,0> 



1 M 1 n -4 OS 



-2t 



2t 



-9.51 
-7.51 
-7.51 
-5.51 



N c =3 m s = i r(r+l)=f s(s+l) = 



_ 3 

4 



f,-l,0,0, -4) 



1 i 

2 ' 2 ' 

1 1 15 

2 ' 2 ' 4 



|,-i,-2i,0,0> 
, f ,i,2i,0,0) 
,f,l,0,4,0) 



2t + t/ 
[/ 
17 

-2f + [/ 



-0.51 
-2.51 
-2.51 
-4.51 
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TABLE V. Eigenkets and eigenvalues for N e =3, spin-down states. 



Number 


Eigenstate 


Energy eigenvalue Value for U/t=5, /i=U/2, and h=0.01t 


N c =3 m s =-i r(r+l) = | s(s+l) = § 



1 

2 ' ' 
1 

2 ' " 
1 

2 ' ' 
1 

2 ' ' 
1 

2 ' ' 



1 

2 ' " 



1 

2 ' " 



1 

2 ' " 



1 

2 ' " 



1 

2 ' ' 
1 

2 ' " 
1 

2 ' " 
1 

2 ' " 



1 

2 ' ' 
1 

2 ' ' 
1 

2 ' ' 



f,-l,-iV3,0,0) 
|,-l,-i>/3,0,0> 
1,-1,1^,0,0) 
1,-1,^,0,0) 



2 




2 




2 + 


a/16* 2 


-4tU- 


-C/ 2 


2 ~ 




2 




U 


x/l6t 2 


-4t IH 


-C7 2 


2 




2 




£ J. 


\/l6* 2 


-4t [7- 


-t7 2 



},-»,», 0,0) 
},-»,», 0,0) 



— + 



},-t,t,0,0> 
},-»,», 0,0) 



- + 



l,t,-t,0,0> 
},<,-», 0,0) 



— + 



!,»,-*, 0,0) 

|,t,-t,0,0> 



— 4- 

9 I 



^/ 32 t 2 + C7 2 +4 V<34 t 4 +3 t 2 [7^ 
2 

y 32 t 2 + [7 2 +4 V<34 t 4 +3 t 2 E7^ 

y 7 32 t 2 +U 2 -4 v / 64t 4 +3t 2 £/ 2 
2 

\/ 32i 2 + i7 2 -4 v / 64i 4 +3t 2 U"^ 
2 

y 7 32 t 2 + U 2 +4 y64 t 4 +3 t 2 (7^ 
2 

y^ 32 i 2 + (7 2 +4 t 4 +3 t 2 (7^ 

2 

y 32i 2 + (7 2 -4 v / 64t 4 +3t 2 £7^ 
2 

y 32i 2 + (7 2 -4 v / 64t 4 +3t 2 £7^ 



1,1,-^,0,0) 
1,1,-^,0,0) 
|,l,t>/3,0,0> 
|,l,t>/3,0,0> 



t7 


Viet 2 


+4t[/+[/ 2 


2 




2 


2 + 


Vwt 2 


+4tU+U 2 


2 T 




2 


[7 


V 16 t 2 


+4t[/+[/ 2 


2 




2 


f + 




+4t[/+[/ 2 






2 



-7.28129 
-2.69871 
-7.28129 
-2.69871 

-10.0929 

0.112923 

-6.55849 

-3.42151 

-10.0929 

0.112923 

-6.55849 

-3.42151 
-8.89512 
-1.08488 
-8.89512 
-1.08488 



N =3 m s =-i r(r+l) = § s(s+l) = 



1 

2 ' " 
1 

2 ' " 
1 

2 ' " 
1 

2 ' ' 



f ,-1,0,0,4) 



f ,-i,-2t,0,0> 
X .*.2»,0,0> 
f ,1,0,-4,0) 



-2t 



2t 



-9.49 
-7.49 
-7.49 
-5.49 



N c =3 m s 



\ r(r+l) = ^ s(s+l) = f 



l 

2 ' ' 
1 

2 ' ' 
1 

2 ' ' 
1 

2 ' ' 



3 
' 4 ' 

3 
' 4 ' 

3 
' 4 ' 

3 
' 4 ' 



-1,0,0, -4) 
-i,-2z,0,0) 
i,2i,0,0) 
1,0,4,0) 



2t + U 
U 

u 

-2t + U 



-0.49 
-2.49 
-2.49 
-4.49 



Ne=3 m a =-| r(r+l) = § s(s+l) = ^ 



l 

2 ' ' 
1 

2 ' ' 
1 

2 ' ' 
1 

2 ' ' 



-1,0,0,4) 
-i,-2i,0,0) 
i,2i,0,0) 
1,0,-4,0) 



-2t 



2t 



-9.47 
-7.47 
-7.47 
-5.47 
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TABLE VI. Eigenkets and eigenvalues for N e =4, spin-up states. The abbreviation a is the same as in Tab. II 



Number 


Eigenstate 




Energy eigenvalue 


Value for U/t=5, n=\l/2, and h=0.01t 








m o **/*>_l_l\ O c _l_ 1 \ R 

m s — z r^r-(-lj — u s^s-(-lj — o 




n i 

94 


|0,2,0,6,-1,0,0,4) 




n 
U 


-1U.U4 








Ills — r I^I^J.^ — u o^a^J.^ — z 




JO 


|U, 1, U, Z, — 1, — A I U, 




I J 
U 


^ 09 


96 


10 10 2—12? a/2 0) 




I J 


-5.02 


97 


0, 1, 0, 2, — i, 0, 0,0} 




2 2 


-10.7216 


98 


0, 1, 0, 2, —i, 0, 0,0) 




r/ | \/i6 t 2 +u 2 

2 ' 2 


-4.31844 


99 


|0,l,0,2,i,0,0,0) 




(7 V 16 t 2 + U 2 
2 2 


-10.7216 


100 


|0, l,0,2,i,0,0,0) 




[7 | y/l6t*+V 3 


-4.31844 


101 


|0,1,0,2,1,0, -4,0} 




2U 2^4»t 2 + U 2 cos(f ) 
3 3 


-1.67211 


102 


10,1,0,2,1,0, -4,0) 




2(7 2^/48t 2 + (7 2 cos(2-^) 
3 3 


-11.5336 


103 


1 2 1 —4 0) 

U, -L, ^, -L, U, 




2(7 2V48t 2 + (7 2 cos(2±^i) 
3 3 


-6.85429 






N —4 


m — 1 rtrA-l} — sCs4-1 1 — fi 




104 


In i n fi — i n n d\ 

U, 1, u, U, 1, u, u, 




n 










111 v J- I l 1 \^ -L J ^ & I J- 1 zL 






n i 9 9 1 n n /)\ 

|U, 1, A A — 1, U, U, 4) 




/■[ 


Pi ri9 


106 


|0,l,2,2,-»,-2»,0,0) 




2t + E/ 


-3.02 


107 


|0, 1, 2,2, -»,2i, 0,0) 




-2t+C/ 


-7.02 


108 


|0,l,2,2,i,-2i,0,0) 




-2t+U 


-7.02 


109 


|0,l,2,2,i,2i,0,0) 




2t + U 


-3.02 


110 


|0,1,2,2,1,0,4,0) 




U 


-5.02 
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VII. Eigenkets and eigenvalues for N e =4, states with spin-projection 0. The abbreviation a is the same as in Tab. 

3 

l\5, 



TABLE 

and /3 is defined as (3 — arccos(4£ 2 U/ ^ 



3 +VJ )• 



Number 



Eigenstate 



Energy eigenvalue 



Value for U/t=5, /x=U/2, and h=0.01t 



N c =4 m s =0 r(r+l)=0 s(s+l)=0 



10,0,0,0,-1,0,0,4} 

10,0,0,0,-1,0,0,4} 

10,0,0,0,-1,0,0,4} 

j0, 0,0,0, -i, -2 i,0,0) 

|0, 0,0,0, -j,2z,0,0} 

]0,0,0,0,i,-2i,0,0} 

j0,0,0,0,i,2i,0,0} 

10,0,0,0,1,0,4,0} 

10,0,0,0,1,0,4,0} 

10,0,0,0,1,0,4,0} 



u 



u + 
u + 



2^/l6t 2 +U 2 cos(-|) 
2 ^/l6t 2 +U 2 cos(i^) 

vf 

2 v / 16t 2 + (7 2 cos( 2-±£ ) 

-2t + U 
2t + U 
2t + U 



U + 



-2t+U 

2^16 t 2 +U 2 cos(|) 
71 ~ 



U ■ 

u 



2 v / 16t 2 +(7 2 cos(; 



2 v / 16t 2 + (7 2 cos(2-t£) 

75 

P = arccos i±JL — r 



-11.8443 
0.844289 

-4. 

-7. 

-3. 

-3. 

-7. 
1.84429 
-10.8443 

-6. 



N c =4 m s =0 r(r+l)=0 s(s+l)=2 



121 


10,0,0,2, 


-1, -2 i a/2, 0,0) 


122 


10,0,0,2, 


-1, 2 i a/2, 0,0) 


123 


|0,0,0,2, 


-i, 0,0,0} 


124 


10,0,0,2, 


-i, 0,0,0} 


125 


10,0,0,2, 


i, 0,0,0} 


126 


10,0,0,2, 


i, 0,0,0} 


127 


10,0,0,2, 


1,0, -4,0} 


128 


10,0,0,2, 


1,0, -4,0} 


129 


10,0,0,2, 


1,0, -4,0} 



u 
u 



U _ V / 16i 2 +C/ 2 

2 2 

V , y / 16i 2 +C7 2 

2 2 

(7 \J\6t 2 +U 2 



u , a/is^+i/ 2 

2 ' 2 
2£7 , 2 y / 48i 2 + Cf 2 cos(f ) 

3 + 3_ 

2U _ 2 ^48 t 2 + U 2 cos (^2.) 

3 3_ 

2(7 2 y/ 48 t 2 + U 2 cos(2±£) 



-5. 

-5. 
-10.7016 
-4.29844 
-10.7016 
-4.29844 
-1.65211 
-11.5136 
-6.83429 



N e =4 m s =0 r(r+l)=0 s(s+l)=6 



0,0,0,6,-1,0,0,4} 







-10. 



131 


10,0,2,0,-1,0,0, -4} 


132 


10,0,2,0,-1,0,0, -4} 


133 


10,0,2,0,-1,0,0, -4} 


134 


|0,0, 2,0, -i, 0,0,0} 


135 


|0,0, 2,0, -i, 0,0,0} 


136 


|0, 0,2,0,i,0,0,0} 


137 


|0, 0,2,0,i,0,0,0} 


138 


|0, 0,2, 0,1, -2 i a/2, 0,0) 



N c =4 m s =0 r(r+l)=2 s(s+l)=0 



2 y/48t 2 + U 2 cos(^i) 



4(7 
3 



4(7 , 2 y / 48t 2 + a 2 cosjl^) 

3 + 3 

4J7 , 2 V48i 2 + LT 2 cos(^l) 
3 + 3 

3 U _ \Zl6i 2 + (7 2 

2 2 
3J7 , x/l6 t 2 + U 2 

2 ' 2 
3 1/ _ y / 16i 2 -H7 2 

2 2 
3(7 v /l6i 2 + (7 2 



■ + 



u 



-8.34789 
1.5136 
-3.16571 
-5.70156 
0.701562 
-5.70156 
0.701562 
-5. 
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139 


10,0,2,0,1,2^,0,0} 


U 


-5. 


N c =4 m s =0 r(r+l)=2 s(s+l)=2 


140 


10,0,2,2,-1,0,0,4} 


U 


-5. 


141 


|0, 0,2,2,-1,-2*, 0,0) 


2t + U 


-3. 


142 


|0,0,2,2,-i,2i,0,0) 


-2t + U 


-7. 


143 


j0,0,2,2,i,-2i,0,0) 


-2t + U 


-7. 


144 


|0,0,2,2, i, 2 i,0,0) 


2t + U 


-3. 


145 


10,0,2,2,1,0,4,0) 


U 


-5. 


N c =4 m s =0 r(r+l)=6 s(s+l)=0 


146 


10,0,6,0,1,0,4,0} 


2U 





TABLE VIII. Eigenkets and eigenvalues for N c =4, spin-down states. The abbreviation a is 


the same as in Tab. III. 


Number 


Eigenstate 


Energy eigenvalue Value for U/t= 


=5, j«=U/2, and h=0.01t 


N c =4 m s = -l r(r+l)=0 s(s+l)=2 


147 


|0,-1, 0,2,-1,-21^2,0,0) U 


-4.98 


148 


|0,-1, 0,2,-1,21^2,0,0) 


U 


-4.98 


149 


10,-1,0,2, -i, 0,0,0} 


U a/ 16 t 2 + U 2 
2 2 


-10.6816 


150 


jO, -l,0,2,-i,0,0,0) 


u . V 16 * 2 +f 2 

2 ' 2 


-4.27844 


151 


10,-1,0,2, i,0, 0,0) 


u \fint 2 +u 2 

2 2 


-10.6816 


152 


10,-1,0,2, i,0, 0,0) 


u | V 16 ' 2 + c/2 

2 ' 2 


-4.27844 


153 


10,-1,0,2,1,0,-4,0) 


2!7 2^48t 2 + U 2 cos(f ) 
3 + 3 


-1.63211 


154 


10,-1,0,2,1,0,-4,0} 


2U 2y/48t 2 + U 2 cos(2-^) 
3 3 


-11.4936 


155 


10,-1,0,2,1,0,-4,0) 


2U 2^4it 2 + U 2 cos(^2) 
3 3 


-6.81429 


N c =4 m s = -l r(r+l)=0 s(s+l)=6 


156 


10,-1,0,6,-1,0,0,4} 





-9.98 


N c =4 m s = -l r(r+l)=2 s(s+l)=2 


157 


10,-1,2,2,-1,0,0,4} 


u 


-4.98 


158 


|0,-l,2,2,-i,-2i,0,0) 


2t+U 


-2.98 


159 


|0,-l,2,2,-i,2i,0,0) 


-2t + U 


-6.98 


160 


|0, — 1, 2,2,4,-21,0,0) 


-2t + U 


-6.98 


161 


|0 I -1,2,2,»,2»,0,0> 


2t + U 


-2.98 


162 


10,-1,2,2,1,0,4,0) 


U 


-4.98 


N c =4 m s = -2 r(r+l)=0 s(s+l)=6 


163 


10,-2,0,6,-1,0,0,4) 





-9.96 
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TABLE IX. Eigenkets and eigenvalues for N e =5, spin-up states. 



Number 


Eigenstate 




Energy eigenvalue Value for U/t=5, 


A*=U/2, and h=0.01t 




N e 


=5 


m s = | r(r+l) = | s(s+l) = f 




164 


1 3 3 15 _i a (| il 




2t + U 


-5.53 


165 


lif, |,f ,-i,2t,0,0> 




U 


-7.53 


166 


lif, |,f ,t,-2t,0,0> 




u 


-7.53 


167 


|I 3 3 15 1 Q _A Q\ 

1 2 ' 2 ' 4 ' 4 ' ' ' ' / 




-2t + U 


-9.53 






=5 


m s = i r(r+l) = | s(s+l) = | 




168 


|i,l,ii-l,-iV^,0,0) 




/ 

3U ^J16t 2 +AtU+U 2 
2 2 


-8.91512 


169 


||,l,i|,-l,-zV3,0,0) 




^ jj a/ 16 f 2 +4 £ U-\-U 2 
2 "r 2 


-1.10488 


1 70 
1 1 U 


I 2' 2 ' 4' 4' — X > 1 V O, U, U; 




3(7 V /l6t2 + 4tc/ + ,;/2 
2 2 


-o.y 10 iz 


171 


1 1 3 1-1,1^,0,0) 
'2' 2' 4 ' 4' ' v ' ' / 




3(7 , V 16 * z +4tU+U z 
2 2 


-1.10488 


172 


i 1 1 3 3 i i Q) 

1 2 ' 2 ' 4 ' 4 ' ' ' ' / 




3(7 -\/32t2 + i72 +4A y 64 j4 + 3 t 2 V 2 
2 2 


-10.1129 


173 


|i, !,§,§,-;,-*, 0,0) 




3(7 | yJz2t 2 + U 2 +4\/mt i +3t 2 U 2 
2 + 2 


0.0929233 


174 


11,1 33 i, i,0,0) 




3!7 i v /32t 2 + (7 2 -4 v / 64t4+3i 2 !7 2 


-6.57849 


1 2 ' 2 ' 4 ' 4 ' ' J "J 




2 2 


175 
176 

1 77 
LI I 


|1 1 3 3 i _i V 

1 2 ' 2 ' 4 ' 4 ' ' J "J 

(2' 2' 4' 4'^'^'^'^^ 

1 1 1 3 3 ■ ■ n n\ 
1 2 ' 2 ' 4 ' 4 ' ' ' ' / 




3!7 ! -\/32t 2 + i7 2 -4 v / 64i 4 +3i 2 U 2 
2 ' 2 

3(7 yj32t 2 + U 2 +i v / 64i 4 +3t 2 (7 2 
2 2 

3(7 \/32t 2 + i7 2 +4 v / 64i 4 +3t 2 (7 2 

-I- — 

2 ' 2 


-3.44151 
-10.1129 

n nnonoQQ 
U.Uyzyzoo 


178 


II I 3 3 ■ " Q Q\ 

1 2 ' 2 ' 4 ' 4 ' ' ' ' / 




3 "\Z 32 t 2 + (7 2 — 4 a/ 64 £ 4 +3 i 2 (7 2 
2 2 


-6.57849 


179 


ll I 2 3 • • Q Q\ 

1 2 ' 2 ' 4 ' 4 ' ' > ' / 




3J7 | -\/32t 2 + i7 2 -4 y / 64i 4 +3i 2 (7 2 
2 "•" 2 


-3.44151 


180 


||,i|,|,i,-»v^,o,o> 




3(7 V /l6 * 2 - 4t!7 + c/2 
2 2 


-7.30129 


181 


|i,l,i|,l,-iV3,0,0) 




3(7 | V 16 ' 2 -*tU+U 2 
2 ' 2 


-2.71871 


182 


|±, i, if,i,^,o,o> 




3(7 y/l6t 2 -4tU+U 2 
2 2 


-7.30129 


183 


|I, I, §,§,1,^,0,0) 




3(7 , \/l6t 2 -4tU+U 2 
2 ' 2 


-2.71871 






=5 


m s = i r(r+l) = § s(s+l) = f 




184 


1 2 ' 2 ' 4 ' 4 ' > ' ' / 




2t + U 


-5.51 


185 


1 1 I 3 15 _ • 2 • q q\ 
1 2 ' 2 ' 4 ' 4 ' '» u ' u / 




U 


-7.51 


186 


|i,l,|,f ,j,-2i,0,0) 




u 


-7.51 


187 


|I 1 2 15 1 n _4 Q\ 

2 ' 2 ' 4 ' 4 ' ' ' ' / 




-2t + U 


-9.51 




No 


=5 


m s = l r(r+l) = f s(s+l) = § 




188 


ll 1 M 2 _i -4) 

1 2 ' 2 ' 4 ' 4 ' ' ' ' / 




-2i + 2(7 


-4.51 


189 


|i,l,f > |,-t,2t,0,0> 




2U 


-2.51 


190 


|i,l,f ,|,j,-2i,0,0) 




2U 


-2.51 


191 


II I 15 1 J n 4 n\ 

1 2 ' 2 ' 4 ' 4 ' ' ' ' / 




2t + 2U 


-0.51 
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TABLE X. Eigenkets and eigenvalues for N e =5, spin-down states. 



Number Eigenstate Energy eigenvalue Value for U/t=5, /tt=U/2, and h=0.01t 











=5 m 


-_i r(r+l) = | s(s+l) = | 




192 


i 

2 ' 


1 

2 


|,|,-i,-iv^,o,o> 




3(7 ^16*2+4* C/+C/ 2 
2 2 


-8.89512 


193 


1 

2 ' 


1 

2 


|,|,-l,-i>/3,0,0> 




3 U i y 10 t u + u ^ 
2 + 2 


-1.08488 


194 


1 

2 ' 


1 

— 2 


I, i V3,0,0> 




3U y/l6t 2 +4tU+U 2 
2 2 


-8.89512 




1 

2 ' 


1 

2 


33 i i ,/i n n\ 
4 1 4 ) x , 1 V o, u, u; 




3(7 | v /l6 * 2 + 4i C/ + Lr2 
2 + 2 


1 nsztss 

- l.UO^OO 


1 Qfi 


1 

2 ' 


1 

2 






3r/ -\/32i 2 + !7 2 +4 v / 64t 4 +3t 2 U 2 
2 2 


1 f) DQ9Q 


197 


1 


1 


33 i i 01 




3(7 | \f&t 2 +U 2 +4,^64,t*+3t 2 U 2 


0.112923 


2 ' 


2 


4545 7 , , / 




2 "•" 2 


198 


1 


1 


?, #, — i, — i, 0, 0> 




3(7 -\/ / 32t 2 + !7 2 -4 v / 64t 4 +3i 2 !7 2 


-6.55849 


2 ' 


2 


4 ' 4 ' ' 5 5 w / 




2 2 


199 


1 


1 


33 2 i 0) 




3(7 ! ^32 t 2 + (7 2 -4 ^/64t 4 +3 i 2 !7 2 


-3.42151 


2 ' 


2 


454) fc 7 fc 7 U 7 U / 




2 ' 2 


200 


1 \ 

2 ' 


_ 2 


f , §,i,i, 0,0) 




3(/ ■\/ 32i2 + !72 + 4 \/ 64t4 + 3 * 2 c/2 
2 2 


-10.0929 


201 


1 

2 ' 


1 

2 


f,f,M,0,0) 




~ Tr 1 / 30 + 2 -X-TT 2 -X-A * / P.A rr2 

3 (7 A/ OiJ t V T^'i y U*± L Tu t- 1^ 

2 2 


0.112923 


202 


1 

2 ' 


1 

2 


f,f,M,0,0) 




3 y y 32t 2 +{7 2 -4y / 64t 4 +3t 2 {7 2 
2 2 


-6.55849 


203 


1 

2 ' 


1 

2 


!,!,<, t, o,o> 




3(7 ! y 32t 2 +E7 2 -4-y/64t 4 +3t 2 E7 2 
2 "r 2 


-3.42151 


204 


1 

2 ' 


1 

2 


f, f,l,-iV3,0,0) 




3(7 v / 16i 2 -4t(7+(7 2 
2 2 


-7.28129 


205 


1 

2 ' 


1 

2 


|,|,l,-tV3,0,0> 




3(7 V /l6 ' 2 - 4t,;/ + Lr2 
2 "r 2 


-2.69871 


206 


1 

2 ' 


1 

2 


|, |, l,i V3,0,0) 




3(7 v /l6 * 2 - 4tc/ + Lr2 
2 2 


-7.28129 


207 


1 


__L 


|,|,l,i^,0,0) 




3(7 | ^/lfit 2 -4tU+U 2 


-2.69871 










=5 m s 


=-i r(r+l) = | s(s+l) = f 




208 


1 

2 ' 


1 

2 


|, f ,-1,0,0,4) 




2t + [/ 


-5.49 


209 


1 

2 ' 


1 

2 


f,i£ ,-t,2t,0,0> 




t/ 


-7.49 


210 


1 

2 ' 


1 

2 


|,f ,i,-2i,0,0) 




(7 


-7.49 


211 


1 

2 ' 


__L 


3 15 1 -4 0) 




— 2i + C/ 


-9.49 










=5 m s 


=-± r(r+l)=f s(s+l) = | 




212 


1 

2 ' 


1 

2 


f ,§,-1,0,0,-4) 




-2i + 2(7 


-4.49 


213 


1 

2 ' 


1 

2 


f ,| ) -i,2i,0,0> 




2(7 


-2.49 


214 


1 

2 ' 


1 

2 


f ,|,i,-2i,0,0) 




2(7 


-2.49 


215 


1 

2 ' 


1 

2 


f ,1,1,0,4, 0) 




2t + 2U 


-0.49 








No 


=5 m s 


=-§ r(r+l) = | s(s+l) = f 




216 


1 

2 ' 


3 
2 


f,f ,-1,0,0,4) 




2i + (7 


-5.47 


217 


1 

2 ' 


3 
2 


|,f ,-i,2t,0,0> 




t/ 


-7.47 


218 


1 

2 ' 


3 
2 


f,f ,1,-21,0,0) 




(7 


-7.47 


219 


1 

2 ' 


3 
2 


2,15,1,0,-4, 0) 




— 2i + C/ 


-9.47 
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TABLE XI. Eigenkets and eigenvalues for N e =6. The abbreviation a is the same as in Tab. 



III. 



Number 


Ei 


genstate 


Energy eigenvalue Value for U/t=5, 


H=U/2, and h=0.01t 


N c =6 m s = l r(r+l)=2 s(s+l)=2 


220 


1 


1,2,2,-1,0,0, -4) 


2U 


-5.02 


221 


1 


1, 2, 2, -i, -2 i,0,0) 


~2t + 2U 


-7.02 


222 


|1, 


1, 2, 2, -i, 2 i,0,0) 


2t + 2U 


-3.02 


223 


|1, 


1,2,2, i,-2i, 0,0) 


2t + 2U 


-3.02 


224 


|1, 


l,2,2,i,2i,0,0) 


-2t + 2U 


-7.02 


225 


|1, 


1,2,2,1,0,-4,0) 


2U 


-5.02 


N c =6 m s =0 r(r+l)=2 s(s+l)=0 


226 


1 


0,2,0,-l,-2iV2,0,0) 


2U 


-5. 


227 


1 


0,2,0,-1,21^,0,0) 


2U 


-5. 


228 


1 


0, 2,0, -i, 0,0,0) 


2 2 


-5.70156 


229 


1 


0,2,0,-i,0,0,0) 


2 ' 2 


0.701562 


230 


1 


0,2,0,i,0,0,0) 


5(7 V 16 * 5 ^ 2 
2 2 


-5.70156 


231 


1 


0,2,0,i,0,0,0) 


5(7 | \/l6t 2 +C/ 2 
2 ' 2 


0.701562 


232 


1 


0,2,0,1,0,4,0) 


7 [7 2- v /48t 2 +(/ 2 cos(f ) 
3 3 


-8.34789 


233 


1 


0,2,0,1,0,4,0) 


7U 2 v /48t 2 + C/ 2 cos(2-jp) 
3 + 3 


1.5136 


234 


1 


0,2,0,1,0,4,0) 


7U 2 1 /48i 2 + U 2 cos (2±H) 
3 3 


-3.16571 


N c =6 m s =0 r(r+l)=2 s(s+l)=2 


235 


1 


0,2,2,-1,0,0, -4) 


2U 


-5. 


236 


1 


0, 2, 2, -i, -2 i,0,0) 


~2t + 2U 


-7. 


237 


1 


0, 2, 2, -i, 2 i,0,0) 


2t + 2U 


-3. 


238 


1 


0,2,2,i,-2i,0,0) 


2t + 2U 


-3. 


239 


1 


0,2,2,i,2i,0,0) 


-2t + 2U 


-7. 


240 


1 


0,2,2,1,0,-4,0) 


2U 


-5. 


N c =6 m s =0 r(r+l)=6 s(s+l)=0 


241 


1 


0,6,0,-1,0,0, -4) 


3U 





N c =6 m s = -l r(r+l)=2 s(s+l)=2 


242 


1 


-1,2,2,-1,0,0, -4) 


2U 


-4.98 


243 


1 


-1, 2, 2, -i, -2 z,0,0) 


-2t + 2U 


-6.98 


244 


|1, 


-1,2,2, -i,2i, 0,0) 


2t + 2U 


-2.98 


245 


1 


-l,2,2,i,-2i,0,0) 


2t + 2U 


-2.98 


246 


1 


-1, 2, 2, i, 2 i,0,0) 


~2t + 2U 


-6.98 


247 


1 


-1,2,2,1,0,-4,0) 


2U 


-4.98 
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TABLE XII. Eigenkets and eigenvalues for N e =7. 



Number 


Eigenstate 


Energy eigenvalue 


Value for U/t=5, /U=U/2, and h=0.01t 




N e 


=7 m s = i r(r+l) = J| s(s+l) = 


3 
4 


248 


111 M 2 -i o -4) 


2t + 3(7 


-0.51 


249 


|§,±,f ,|,-t,-2i,0,0> 


3(7 


-2.51 


250 


|§,±,±5 ,|,t,2i,0,0> 


3(7 


-2.51 


251 


|3 I 15 3 x q 4 

1 2 ' 2 ' 4 ' 4 ' ' ' ' / 


-2t + 3U 


-4.51 




N c = 


=7 m s =-i r(r+l)=f s(s+l) = 


_ 3 

4 


252 


|3 -I 3 _ 1 q q _ 4 \ 


2t + 3U 


-0.49 


253 


1 2 ' 2 ' 4 ' 4 ' ^' 2 2, 0, 0} 


3(7 


-2.49 


254 


||,-I,^,|,t,2i,0,0> 


3(7 


-2.49 


255 


1 2 ' 2 ' 4 ' 4 ' ' ' ' / 


-2t + 3U 


-4.49 



TABLE XIII. Eigenkets and eigenvalues for N c =8. 



Number 


Eigenstate 


Energy eigenvalue Value for U/t=5, /i=U/2, and h=0.01t 


N c =8 m s =0 r(r+l)=6 s(s+l)=0 


256 


12,0,6,0,1,0,4, 0) 


4(7 
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. . . . 4^ 



4^ 




Figure 1 



U=5t, <N e >=4, h=0.01t 



n=5, r=1/2, s=3/2, u=1 



n=3, r=1/2, s=3/2, u=-1 



n=4, r=0, s=2, u=-1 



n={3,5}, r=1/2, s=1/2, u={i,-i} 



n=4, r=0, s=1 , u={i,-i} 



n=4, r=0, s=0, u=1 



n=4, r=0, s=1 , u=-1 



n=4, r=0, s=0, u=-1 



Figure 2 
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n=3, r=1 /2, s=1/2, u=1 , k 1 ={3 ,/2 i,-3 1 ' 2 



n=4, r=0, s=1 , u={i,-i} 
. n=4, r=0, s=0, u=1 



n=3, r=1/2, s=3/2, u=-1 



n=4, r=0, s=0, u=-1 



n n=4 


- — t 


n=3 


- —i 


n-3 n=2 


- / 





Figure 5 
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46, 50, 70, 74 : 
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-4.75 
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1 1 1 , 22 
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Figure 8 
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